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Abstract 

A measurement of inclusive jet cross-sections in deep-inelastic ep scattering at HERA 
is presented based on data with an integrated luminosity of 21.1 pb _1 . The measurement 
is performed for photon virtualities Q 2 between 5 and 100 GeV 2 , differentially in Q 2 , in 
the jet transverse energy Et, in E T /Q 2 and in the pseudorapidity 77i a b- With the renor- 
malization scale [Mr = Et, perturbative QCD calculations in next-to-leading order (NLO) 
give a good description of the data in most of the phase space. Significant discrepancies are 
observed only for jets in the proton beam direction with Et below 20 GeV and Q 2 below 
20 GeV 2 . This corresponds to the region in which NLO corrections are largest and further 
improvement of the calculations is thus of particular interest. 
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1 Introduction 



HERA offers excellent possibilities to test predictions of quantum chromodynamics (QCD), the 
theory of the strong interaction, in deep-inelastic electron^ proton scattering (DIS). This theory 
has been successfully tested to a very high level of precision in the measurements of the proton 
structure function F 2 [0],@], in which the photon virtuality Q 2 gives rise to a hard scale. QCD 
also predicts the production of partons with large transverse momenta, which fragment into 
hadronic jets with similar four-momenta. Jet observables therefore give direct access to the un- 
derlying parton dynamics. In addition to ■JQ*, the transverse energies E T of the resulting jets, 
measured in an appropriate frame of reference, provide a natural hard scale for the description 
of the interaction within perturbative QCD. 

Jets have been studied extensively at HERA and at other colliders. Recent measurements 
of (multi-)jet cross-sections at HERA have shown that at sufficiently large values of the 
photon virtuality Q 2 ^150 GeV 2 , calculations implementing matrix elements to 0(a 2 ) in the 
strong coupling constant, i.e. next-to-leading order (NLO), are in excellent agreement with the 
data. These highly successful calculations are performed in the conventional collinear approx- 
imation referred to as DGLAP [0]. Since in this kinematic regime of large Q 2 , the theoretical 
uncertainties are small, the measured cross-sections have been used to extract a s [|3|,|5p and the 
gluon density of the proton [|3|]. For some observables these data can be described even at lower 
values of Q 2 , although the theoretical uncertainties in the predictions and the QCD corrections 
between leading order (LO) and NLO predictions increase in this kinematic region. 

Problems have also been encountered in the theoretical description of jets in DIS. Measure- 
ments of highly energetic jets at large pseudorapiditiesf] ?7i ab with transverse energies squared of 
the order of Q 2 (so-called forward jets [§) made by the HI [§] and ZEUS [0|T|] collabora- 
tions show large discrepancies between data and NLO calculations at low values of Bjorken-x. 
This may be due to a breakdown of the DGLAP approximation and the onset of BFKL [ [T2TJ 
effects where, in contrast to DGLAP, the emissions of partons are not ordered in transverse 
momentum. Alternatively, a good description of the data can be obtained [ |l3l ] by considering 
the partonic structure of the virtual photon, which is expected to be important for E\ > Q 2 
and Q 2 not too large. It has been pointed out Jl4|], however, that there is some correspondence 
between higher order effects and resolved virtual photon contributions. Further discrepancies 
between NLO calculations and data were revealed by HI measurements of DIS dijet production 
for Q 2 < 11 GeV 2 |J5p and for Q 2 < 10 GeV 2 with invariant dijet masses M sj < 25 GeV or 
mean dijet E T < 20 GeV [|3p. In these regions of phase space the NLO corrections to the LO 
cross-sections are large, such that next-to-next-to-leading order (NNLO) contributions are also 
likely to be important. 

While the HI measurements of dijet [|3p and three-jet production [Q] mentioned above cover 
low and high Q 2 values, inclusive jet production in DIS has only been measured previously at 
high values of Q 2 above 150 GeV 2 [J3p. This paper presents precision measurements of inclusive 

1 In the following, the generic name "electron" will be used for the beam and scattered lepton. 

2 The pseudorapidity rji^b is defined as ?yi a b = — lntan(0/2), with the polar angle 6 being measured with 
respect to the positive z-axis, which is given by the proton beam (or forward) direction. 
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jet cross-sections in the Breit frame[] at low Q 2 values (5 < Q 2 < 100 GeV 2 ). Here, "inclu- 
sive" means that in every selected event every jet that passes the experimental cuts contributes 
to the measured cross-section. The advantages of investigating QCD with studies of inclusive 
jets compared with multi-jets are two-fold. First, the accessible phase space is extended. Sec- 
ond, phase space regions in which fixed order calculations are infrared- sensitive are naturally 
avoided. In dijet analyses such regions are usually suppressed by imposing asymmetric cuts on 
the transverse jet energies fifflft. 



2 Experimental apparatus 

A detailed description of the HI apparatus can be found elsewhere [|17|,[l8[]. Here, only the 
components of the detector which are relevant for the measurement of low Q 2 jet cross-sections 
are introduced. 

The scattered electron is measured with a calorimeter made of lead and scintillating fibres 



(SpaCal) fll9Q. The SpaCal covers polar angles from 153° to 177.5° for collisions at the nominal 
interaction point. It is divided into an electromagnetic section with 28 radiation lengths and a 
hadronic section and has a total depth of 2 interaction lengths. The electromagnetic section of 



the SpaCal has an energy resolution of a E /E ps 0.07/ A/-E[GeV] © 0.01 for electrons [gOp. In 
addition, a backward planar drift chamber in front of the SpaCal with an angular acceptance 
of 151° < 9 < 177.5° serves to suppress photoproduction background, where a high energy 
hadron fakes an electron. 

The hadronic final state is measured using the SpaCal and the liquid argon (LAr) calorime- 
ter pi) ] together with the tracking chambers. The LAr calorimeter provides full azimuthal 
coverage over the polar angle range 4° < 9 < 154° with a depth ranging between 4.5 and 8 in- 
teraction lengths, depending on the polar angle. Test beam measurements of the LAr calorimeter 



modules showed an energy resolution of <Je/E « 0.50/ y £7 [GeV] © 0.02 for charged pions af- 



ter software energy reweighting [ |2"2| ] and of a E /E « 0. 12/ E[GeV] ©0.01 for electrons [|2~3|]. 

The calorimeters are surrounded by a superconducting solenoid which provides a uniform 
magnetic field of 1.15 T parallel to the beam axis in the region of the central tracking detec- 
tors. Charged particles are measured in this central tracking area which consists of drift and 
proportional chambers. The drift chambers cover a range in polar angle from 15° to 165°. The 
tracking chambers also serve for the determination of the event vertex. 



3 QCD calculations and Monte Carlo models 



For the NLO calculation of jet observables at the parton level the DISENT computer pro- 



gram \ yM was used. Several other programs [g5J-g7J are known to give comparable results [|28|]. 
All of these programs calculate the direct photon induced contributions to the cross-section. 

3 The Breit frame is defined by 2xp + q = 0, where x is the Bjorken scaling variable, and p and q are the proton 
and the virtual photon momenta, respectively. 
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Only the JetViP program [ [251 ] provides the additional possibility to calculate a cross-section 
consisting of both direct and resolved photon contributions. While implementing the concept of 
the resolved hadronic substructure of the photon is straightforward for photoproduction, in the 
DIS case conceptual difficulties are encountered [ ^9[ ] which lead to ambiguous results. There- 
fore we decided not to use JetViP for comparisons with the present data. 

In the DISENT program the square of the renormalization scale, fi R , can be set to any linear 
combination of the two hard scales in the event, E\ and Q 2 . Here, Ex is the transverse energy 
of the jet in the Breit frame in the case of events in which only one jet satisfies the selection 
criteria. In the case of two or more selected jets, E T was set to the mean transverse energy of the 
two hardest jets. For most of the comparisons with data fx R = E T was chosen, since E\ 
in almost all of the phase space considered. However, the effects of choosing fx R = 
also discussed in section |5j. The factorization scale was always set to \[Q 2 . In the kinematic 
region considered in this analysis the cross-section predictions are stable within a few percent 
even for large variations of the factorization scale [ ]30| ] . The parton density functions (PDF) of 
the proton were taken from the CTEQ5M (CTEQ5L) parameterization [^TJ for the calculation 
of the NLO (LO) cross-sections. For NLO, the corresponding value of a s (M z ) is 0.118. The 
number of active flavours was chosen to be n/ = 5 unless fx 2 R was less than 25 GeV 2 in which 
case nj = 4 was used in order to simulate the threshold for the onset of beauty production via 
the boson gluon fusion process. 

In contrast to fixed order QCD calculations, which can only predict the partonic final state 
of an event, the implementation of phenomenological QCD models in Monte Carlo (MC) gen- 
erators allows the details of the hadronic final state to be simulated. These generators typically 
implement 0(a s ) matrix elements and account for the effects of higher orders using different 
QCD-inspired mechanisms. LEPTO @ and RAPGAP @ are so-called ME+PS Monte Carlo 
programs which combine the 0(a s ) matrix elements with parton showers which take into ac- 
count the leading logarithms in Q 2 to all orders. As an alternative approach, the ARIADNE 
MC [Ql simulates parton cascades as a chain of independently radiating colour dipoles accord- 
ing to [|35J]. Higher order QED corrections, which can influence the event topology and the size 
of the cross-sections, are implemented in the program HERACLES [j36[] which is interfaced 



to RAPGAP. An interface to LEPTO and ARIADNE is provided by the DJANGO []37|] pro- 
gram. In the present analysis the MC generators ARIADNE and LEPTO were used to estimate 
the hadronization corrections which were applied to the NLO QCD calculations. RAPGAP and 
the DJANGO-interfaced ARIADNE program ( ' DJANGO/ ARIADNE ' ) were used to correct the 
data for detector and higher order QED effects (see section [T2j). For all MC generators used in 
the analysis, the hadronization of the final partonic system as well as the particle decays were 
modelled with the Lund colour string model [p8|] as implemented in JETSET 



4 Measurement 

4.1 Event and jet selection 

The analysed data were collected with the HI detector in the years 1996 and 1997. In this 
running period HERA collided 27.5 GeV positrons with 820 GeV protons. The integrated 
luminosity as measured using the bremsstrahlung process ep — > ep^y amounts to 21.1 pb -1 . 
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The events are triggered by demanding a localized energy deposition in the SpaCal and 
loose track requirements. The trigger efficiencies are close to 100 %. DIS events are selected 
by identifying the scattered electron in the SpaCal. More precisely, the electron is defined as 
the highest energy cluster in the SpaCal with an energy of at least 10 GeV. It must be detected 
at a polar angle 9 > 156° to guarantee the full reconstruction of the particle within the SpaCal 
acceptance. 

The kinematic variables x, Q 2 and y are determined using electron information only, ac- 
cording to Q 2 = AE e E' e cos 2 (fl e /2), y = l- (E' e /E e ) sin 2 (# e /2) and x = Q 2 /(ys), where E e , 
E' e and 9 e are the electron beam energy, the energy of the scattered electron and the electron 
scattering angle, respectively. The latter two quantities are measured using the SpaCal. The 
center-of-mass energy squared s of the proton-electron scattering is computed from the beam 
energies. 

The kinematic range considered in this analysis is further constrained by the conditions 
5 < Q 2 < 100 GeV 2 and 0.2 < y < 0.6. The latter condition also leads to a reduction 
of photoproduction background and of events in which the incoming electron radiates a high- 
energy photon. 

The hadronic final state objects of an event are reconstructed using information from energy 
deposits in the calorimeters and from tracks in the inner detectors. A cut ~ Pz,j) > 

45 GeV further suppresses radiative events. Here the sum runs over the objects of the hadronic 
final state and the scattered lepton. An additional cut J2j(Ej — P z j) < 65 GeV reduces non-ep 
background. For all events the reconstructed ^-coordinate of the event vertex is required to be 
within ±35 cm of its nominal position which further reduces non-ep background to a negligible 
level. 



Jets are defined using the inclusive k± cluster algorithm p0[] in the Breit frame in which the 



photon and the parton collide head-on. This ensures that the jet transverse energies are closely 
related to those of the partons emerging from the hard scattering. Jets are finally selected by 
requiring their transverse energy to be larger than 5 GeV. In addition, the jets have to be well 
contained within the polar angle acceptance of the LAr calorimeter. Therefore, a cut on the 
pseudorapidity of the jets —1 < ^i ab < 2.8 is applied. 



4.2 Correction procedure 

All data distributions shown in this paper are corrected for the effects of limited detector ac- 
ceptance and resolution and for higher-order QED effects using a bin-to-bin correction method. 
Following studies of migrations made with simulated events, the bin sizes were chosen in order 
to ensure that this method can be applied. In all bins, the stability^ and purity exceed 40%. The 
distributions of kinematic variables and jet observables in the data are sufficiently well described 
by the QCD Monte Carlo models used in the correction procedure. RAPGAP describes the data 
very well for transverse energies above 10 GeV but lies below the data for E T < 10 GeV. The 
DJANGO/ ARIADNE program gives a good description for E T < 10 GeV, but is too high for 
higher values of the transverse energy. 

4 The stability (purity) is defined as the number of jets which are both generated and reconstructed in an analysis 
bin, divided by the total number of jets that are generated (reconstructed) in that bin. 
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The correction factor is defined as the ratio of the number of events from two different MC 
event samples generated for the same phase space. The first sample has no QED corrections 
and is not subjected to the detector simulation. The second sample includes QED corrections 
to the leptonic vertex of the interaction and is subjected to a full simulation of the HI detector 
and reconstruction procedure. The correction factors were calculated from both the RAPGAP 
and the DJANGO/ ARIADNE programs. The final correction was taken to be the average of 
the two results. The corrections factors thus derived are typically between 0.9 and 1.1. Val- 
ues down to 0.7 and up to 1.3 are reached in some exceptional analysis bins, for example for 
Et/Q 2 < 2 and E\jQ 2 > 50 for jets in the forward region 1.5 < 7? lab < 2.8. Although the 
shapes and normalizations of the predicted cross-sections differ significantly between RAGPAP 
and DJANGO/ ARIADNE, the difference between the correction factors derived with the two 
programs is usually as small as 5 to 10%. Half of this difference is assigned as the systematic 
error in the correction procedure. 

In order to compare the data, which are corrected for instrumental and radiative effects, 
to the NLO QCD calculations, hadronization corrections are applied to the latter. They were 
estimated using the predictions of the Monte Carlo programs LEPTO and ARIADNE for the 
ratio of the hadron level to the parton level cross-sections. The average of the two predictions is 
taken as the correction. The effects of hadronization reduce the NLO cross-sections by factors 
that vary from 0.85 to 0.95 with increasing jet transverse energy or Q 2 pTIJ. The uncertainty 
in the hadronization correction is taken to be half the difference between the predictions of the 
two programs. 

The treatment of systematic uncertainties is similar to that in [0,0]. The dominant con- 
tribution to the total systematic error stems from the uncertainty in the hadronic energy scale 
of the LAr calorimeter, which is taken to be 3% for this analysis. This uncertainty leads to 
an uncertainty in the measured cross-sections of 10 to 15%. All other sources of systematic 
uncertainties, such as the method used for the reconstruction of kinematic variables, the deter- 
mination of the energy and polar angle of the scattered electron and the model dependence of 
the correction procedure, are small by comparison. 



5 Results 



In this section the results of the inclusive jet measurement described in section |] are com- 
pared with LO and NLO QCD calculations as explained in section ||. The data are presented 
in figures [j] to 0. In the top parts of the figures, the data, corrected for detector and radiative 
effects, are compared with the results of the LO and NLO QCD calculations without hadroniza- 
tion corrections. The uncertainty in the NLO predictions, which was estimated by varying the 
renormalization scale by a factor of ±2, is also shown. In the bottom part of the figures, 
the consistency of the NLO QCD prediction with the data is studied in more detail by plotting 
the relative difference (QCD-Data)/Data. Here "QCD" denotes the NLO prediction corrected 



for hadronization effects as discussed in section 4.2. Both the uncertainty due to the hadroniza 



tion correction and the effect of the variation of the renormalization scale \i R are indicated. All 
measured cross-sections and their errors are given in tables [II |2| and 0, together with the average 
Bjorken-x and Q 2 values of the measurements, taken from the data. 



8 



Figure |T] shows the inclusive jet cross-section as a function of the transverse jet energy E T 
in different regions of the pseudorapidity r]i ab : in the backward region — 1 < r] lab < 0.5, the 
central region 0.5 < ??i a b < 1.5 and the forward (proton direction) region 1.5 < ?7i a b < 2.8. 
The measured cross-sections, which extend over four orders of magnitude, are compared to the 
QCD calculations using a renormalization scale = Et- 

While there is good agreement between the data and the NLO QCD calculation in the back- 
ward region for all Et values, discrepancies are observed for more forward jets with low Et. 
At the lowest Et, 5 < Et < 20 GeV, for r^b > 1.5, the assumed renormalization scale uncer- 
tainty does not cover the large difference between data and QCD calculation. Figure [I] shows 
that the discrepancies are accompanied by large corrections between LO and NLO calculations. 
The NLO predictions are up to a factor 5 larger than the LO ones. Furthermore, the effects of a 
variation of the renormalization scale, which are indicated by the hatched band, are largest for 
these low E T jets in the forward region. 

In figure the jet E T distribution in the forward region only is presented in five different 
intervals of Q 2 . Whereas for Q 2 > 20 GeV 2 the data can be described by the NLO predictions 
even for transverse energies below 20 GeV, for Q 2 < 20 GeV 2 the NLO calculation falls short 
of the data by a factor of up to two for E T < 20 GeV. As in fig. [j] these discrepancies between 
data and predictions are accompanied by large NLO/LO corrections. 

One can conclude that the perturbative NLO QCD calculations work reasonably well in 
most of the rapidity range, even in the forward region i]\ ah > 1.5, as long as both E T and Q 2 are 
not too small. 

In order to study the interplay of the two possible scales in DIS jet events, figure |3] shows 
the inclusive jet cross-section as a function of the ratio E T /Q 2 for the different regions of 
T^iab- The data are well described by the NLO calculations over the full E T /Q 2 range only 
in the backward region —1 < 77^ < 0.5. For the central and even more for the forward 
pseudorapidities, there is a discrepancy between data and calculation for the medium range of 
this ratio, i.e. 2 < E T /Q 2 < 50. This region is dominated by small values of E\ and Q 2 . For 
large and small values of E T /Q 2 , which are correlated with large values of either E T or Q 2 , the 
NLO calculation is in agreement with the data. As in the case of daj et /dE T the discrepancies 
occur in regions where the NLO/LO corrections are very large (up to a factor of 6 in the forward 
region). 

Switching to = \[Q 2 (figure 0) changes the situation. In this case, large discrepancies 
between data and NLO calculation are observed only for large values of E T /Q 2 > 50 (and 
hence low values of Q 2 ) irrespective of the ??i ab range. This suggests that when Q 2 is much 
smaller than E T , the variable Q 2 provides too soft a scale and is irrelevant for the description of 
the hard scattering process. In contrast to the case fi R = E T , these discrepancies do not occur 
in the analysis intervals in which the NLO/LO corrections are largest. Although the data and 
the predictions are compatible at lower values of E T /Q 2 , it should be noted that the choice of 

— \[Q 2 leads to much larger scale uncertainties than is the case with fin = Et- 

6 Summary 

Inclusive jet cross-sections have been measured for values of Q 2 between 5 and 100 GeV 2 with 
the HI detector at HERA. Jets were selected using the inclusive k ± algorithm in the Breit frame 
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and were required to have a minimum transverse energy of 5 GeV. QCD calculations up to 
second order of the strong coupling constant a s , with E T or ^/Q 2 as renormalization scale, 
were tested against the data. 

With E T as renormalization scale, the data on daj ct /dE T are well described by the NLO 
QCD calculations for jets of all transverse energies in the backward region and in all ?7i ab regions 
for high transverse jet energies E T > 20 GeV, even though the NLO cross-sections can be 
more than 3 times larger than the LO cross-sections. The only region in which the NLO QCD 
predictions fail to describe the data is in the forward region when both the transverse jet energy 
and Q 2 are relatively small. In this region, discrepancies of up to a factor of two are observed, 
not considering the scale uncertainties. The regions in which the data cannot be described by the 
QCD calculations are characterized by large NLO/LO corrections and by a strong dependence 
on the renormalization scale. It is worth noting that the theoretical uncertainties, which are 
dominated by the effects of a variation of the renormalization scale, are usually much larger 
than the experimental errors. 

Similarly, the distributions of E\jQ 2 are described by the QCD calculations using jiR = Et 
everywhere except for the range 2 < E\jQ 2 < 50 for jets in the forward and central regions. 
Again, the discrepancies between data and calculation are accompanied by large NLO/LO cor- 
rections. For the choice fiR = \[Q 2 the agreement between data and predictions is in general 
good. However, a very strong dependence of the QCD calculations on the renormalization scale 
is found which questions whether they can be considered as predictive. For [i R = \[Q 2 and 
E\jQ 2 > 50 the calculations fail to describe the data for all jet pseudorapidities, showing that 
Q 2 is not an appropriate scale choice if Q 2 <C E\. 

The correlation of large NLO/LO corrections and high sensitivity to renormalization scale 
variations with poor agreement between data and QCD predictions in the case of = E T 
strongly suggests that the inclusion of higher order (e.g. NNLO) terms in the QCD calculations 
is necessary in order to describe the data. Since in most analysed intervals the experimental er- 
rors are much smaller than the theoretical uncertainties, it is evident that a deeper understanding 
of low Q 2 DIS jet production will require theoretical progress. 
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Figure 1: Inclusive jet cross-sections da ic J6.E T in different ranges of ?7i ab , integrated over the 
region 5 < Q 2 < 100 GeV 2 and 0.2 < y < 0.6. The data are shown as points with error 
bars which include statistical and systematic errors. In the top part the data are compared to 
DISENT NLO QCD calculations using the CTEQ5M parton distribution functions (solid line) 
and to DISENT LO calculations using CTEQ5L (dashed line). The renormalization scale is set 
to /j, R = E T and no hadronization corrections are applied to these predictions. The hatched 
band around the NLO prediction stems from variations of the renormalization scale by factors 
of ±2. In the bottom part the ratio (QCD-Data)/Data is shown. Here "QCD" denotes NLO QCD 
corrected for hadronization effects. The inner white band shows the hadronization uncertainty. 
The hatched outer area contains in addition the renormalization scale uncertainty added linearly. 
A dotted line at serves to guide the eye. 
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Figure 2: Inclusive jet cross-sections daj et /dE T for the forward region 1.5 < 77i ab < 2.8 in 
different ranges of Q 2 . See the caption of figure [l] for further explanation. 



15 



HI Inclusive Jets 



a 10 
o 

tT 10 

73 



10 



as 

i 2 

15 1 
Q 

6 

U 

o 

w -1 



I 



1 '° < ^lab < 05 

m 




i 

i 

HI Data 
LO QCD 



NLO QCD ^ 

M"R =E T 



hadronization 
uncertainty 



°- 5<T llab <1 - 5 



» i 







; 



— s- 




1.5 < Ti lab < 2.8 



mm 



-t— i 



10 



100 1 



10 



100 1 



10 



100 
E 2 /Q 2 



Figure 3: Inclusive jet cross-sections daj ct /d(E^/Q 2 ) in different ranges of ^i a b. See the caption 
of figure |T] for further explanation. 
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dajet/dEr 


»yiab 




value 


stat. error 


syst. error 


< x > 


< Q ' A > 


range 


interval [GeV] 


[pb/GeV] 


[pb/GeV] 


[pb/GeV] 


[lO" 3 ] 


[GeV 2 ] 




5 - 10 


561 


4 


+33 / -35 






-1.0-0.5 


10-20 


31.3 


0.5 


+2.7/ -2.0 


0.7 


21 




20-40 


0.43 


0.03 


+0.06 / -0.06 








5 - 10 


272 


3 


+21/ -19 






0.5 - 1.5 


10-20 


28.2 


0.5 


+2.3/ -2.1 


0.8 


23 




20-40 


1.40 


0.07 


+0.12/ -0.13 








40-70 


0.017 


0.006 


+0.004/ -0.004 








5 - 10 


238 


3 


+18/ -18 






1.5-2.8 


10-20 


21.4 


0.5 


+1.6/ -1.9 


0.8 


21 




20-40 


0.94 


0.06 


+0.07/ -0.13 








40-70 


0.026 


0.008 


+0.004/ -0.006 







Table 1: Measured jet cross-sections daj et /dE T in different ranges of ^i ab for 5 < Q 2 < 
100 GeV 2 and 0.2 < y < 0.6, together with their statistical and systematic uncertainties. The 
mean values of Bjorken-rr and Q 2 of the accepted events are also given. 
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&<j Jct l&E T for 1.5 < ryi a b < 2.8 


Q 


TP 


value 


stat. error 


syst. error 


< x > 


< Q > 


range [ijev J 


interval [GeV] 


LpD/(jev J 


Lpu/<jev J 


LpD/(je v J 


[10 1 


LGeV J 




5- 10 


104 


3 


+9/-11 






5 - 10 


10-20 


7.6 


0.4 


+0.6 / -0.6 


0.2 


7 




20-40 


0.28 


0.03 


+0.05 / -0.02 








40-70 


0.010 


0.005 


+0.003 / -0.002 








5- 10 


68.1 


1.6 


+4.5 / -6.0 






10-20 


10-20 


5.5 


0.2 


+0.4/ -0.4 


0.5 


14 




20-40 


0.25 


0.03 


+0.04/ -0.01 








40-70 


0.005 


0.003 


+0.001 / -0.002 








5- 10 


34.3 


1.4 


+4.1/ -1.6 






20-35 


10-20 


3.7 


0.2 


+0.3 / -0.4 


0.9 


26 




20-40 


0.17 


0.03 


+0.03 / -0.02 








40-70 


0.01 


0.005 


+0.001 / -0.002 








5- 10 


28.0 


0.8 


+1.9/-2.1 






35-70 


10-20 


3.5 


0.2 


+0.2 / -0.4 


1.7 


49 




20-40 


0.16 


0.02 


+0.03/ -0.01 








40-70 


0.005 


0.005 


+0.001 / -0.002 








5- 10 


6.5 


0.4 


+0.5 / -0.8 






70 - 100 


10-20 


1.03 


0.08 


+0.12/ -0.06 


3.4 


82 




20-40 


0.055 


0.015 


+0.008 / -0.006 







Table 2: Measured jet cross-sections daj ct /dET for 0.2 < y < 0.6 in different ranges of Q 2 
for jets in the forward region 1.5 < i]^ < 2.8, together with their statistical and systematic 
uncertainties. The mean values of Bjorken-rr and Q 2 of the accepted events are also given. 
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da Jct /d(E^/Q') 


7?lab 




value 


stat. error 


syst. error 


< x > 


< Q l > 


range 


interval 


[pb] 


[pb] 


[pb] 


[lO" 3 ] 


[GeV 2 ] 




0.5-2 


593 


6 


+32/ -35 


1.2 


37 




2-8 


286 


3 


+17/ -18 


0.4 


14 


-1.0-0.5 


8 - 15 


46.3 


0.9 


+4.7 / -2.7 


0.3 


9 




15-50 


3.3 


0.1 


+0.3 / -0.3 


0.2 


8 




50 - 100 


0.09 


0.01 


+0.02 / -0.01 


0.2 


7 




0.5-2 


280 


4 


+25/ -15 


1.4 


39 




2-8 


149 


2 


+9/ -12 


0.5 


15 


0.5 - 1.5 


8 - 15 


28.9 


0.7 


+2.2/ -2.5 


0.4 


11 




15-50 


3.3 


0.1 


+0.3 / -0.2 


0.3 


11 




50 - 100 


0.22 


0.02 


+0.03 / -0.02 


0.3 


9 




0.5-2 


213 


5 


+18/ -13 


1.4 


38 




2-8 


135 


2 


+10/ -12 


0.5 


14 


1.5-2.8 


8 - 15 


24.5 


0.8 


+2.2/ -2.2 


0.3 


10 




15-50 


2.6 


0.1 


+0.3 / -0.2 


0.3 


9 




50 - 100 


0.16 


0.02 


+0.03 / -0.01 


0.3 


9 



Table 3: Measured jet cross-sections da,] ct /d(E^/Q 2 ) in different ranges of r] lab for 5 < Q 2 < 
100 GeV 2 and 0.2 < y < 0.6, together with their statistical and systematic uncertainties. The 
mean values of Bjorken-x and Q 2 of the accepted events are also given. 
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